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Q^ ■ We consider random systems of equations over the reals, with m equations and ni unknowns 

Pi(i) + Xi(i) = 0, f G R™, i = l,...,m, where the Pi's are non-random polynomials having 
degrees di's (the "signal") and the X^'s (the "noise") are independent real- valued Gaussian cen- 
2 ' tered random polynomial fields defined on R™, with a probability law satisfying some invariance 

properties. 

For each i, Pi and Xi have degree di. 

The problem is the behavior of the number of roots for large m. We prove that under specified 
conditions on the relation signal over noise, which imply that in a certain sense this relation is 
neither too large nor too small, it follows that the quotient between the expected value of the 
VQ ' number of roots of the perturbed system and the expected value corresponding to the centered 

r\l , system (i.e.. Pi identically zero for all i = 1, . . . ,7n), tends to zero geometrically fast as m tends 

f^ ' to infinity. In particular, this means that the behavior of this expected value is governed by the 

I ' ^ noise part. 

>— ' , Keywords: random polynomials; Rice formula; system of random equations 

o 

1. Introduction and main result 

Let /=(/i,...,/„), 

/,(t):= Y. °?^ (* = l,...,m) (1) 

\\3\\<d, 

be a system of m polynomials in m real variables. The notation in (1) is the following: 
i := (ii, . . . ,im) denotes a point in M™, j := (j'l, . . . ,jm) a multi-index of non-negative 
integers, \\j\\ = Y^^=i Jh-, t^ = t^^ ■ ■ ■ P"^ , a^ = aj^ . di is the degree of the polynomial 
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250 D. Armentano and M. Wschebor 

We denote by N^{y) the number of roots of the system of equations 

Mt) = (^ = l,...,m) 

lying in the subset y of M". Also N^' == iV/(R™). 

If we choose at random the coefficients {af }, Nf{V) becomes a random variable. 
Classical results in the case of one polynomial in one variable seem to have started with 
the work of Marc Kac [6] (see the book by Bharucha-Rcid and Sambandham [3]). Here 
we will be interested in systems with m > 1, and more specifically, in large values of m. 
This appears to be of a quite different nature than the case m = 1 and generally speaking, 
little is known on the distribution of the random variable N^ (V) (or N-^) even for simple 
choices of the probability law on the coefficients. 

In 1992, Shub and Smalc [10] computed the expectation of N^ when the coefficients 
are Gaussian centered independent random variables having variances: 

Their result was 

EiNf)^idi---dm)'/^ (3) 

Some extensions of their work, including new results for one polynomial in one variable, 
can be found in [5]. There are also other extensions to multi-homogeneous systems in 
[8], and, partially, to sparse systems in [7] and [9]. A similar question for the number 
of critical points of real-valued polynomial random functions has been considered in a 
recent paper by Dedieu and Malajovich [4]. 

A general formula for E{Nf{V)) when the random functions fi {i — l,...,m) are 
stochastically independent and their law is centered and invariant under the orthogonal 
group on R™ can be found in [1] . This includes the Shub-Smalc formula (3) as a special 
case. Very little in known on higher moments. The only published results of which the 
authors are aware concern asymptotic variances as m — *■ -|-oo (see [1] for non-polynomial 
systems and [11] for the Kostlan-Shub-Smale model). 

The aim of this paper is to remove the hypothesis that the coefficients have zero 
expectation (in some cases, this has been considered for one polynomial in one variable 
in the above-mentioned paper [5]). 

One way to look at this problem is to start with a non-random system 

P,it)^0 (i = l,...,m), (4) 

perturb it with a polynomial noise {Xi(t) :i — 1, . . . , m}, that is, consider 

P^it)+X^{t)^0 {i^l,...,m) 

and ask what one can say about the number of roots of the new system. Of course, to 
obtain results on E{N^'^^) we need a certain number of hypotheses both on the "noise" 
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X and the class of polynomial "signals" P, especially the relation between the size of P 
and the probability distribution of X . 

Roughly speaking, we prove in Theorem 2 that if the relation signal over noise is neither 
too big nor too small, in a sense that will be made precise later on, there exist positive 
constants C,9, < 9 <1 such that 

E{NP+^) < Ce^'EiN^). (5) 

Inequality (5) becomes of interest if the starting non-random system (4) has a large 
number of roots, possibly infinite, and m is large. In this situation, the effect of adding 
polynomial noise is a reduction at a geometric rate of the expected number of roots, as 
compared to the centered case in which all the Pi's are identically zero. 

Notice that in formula (5), E{N^) can be computed by means of a nice formula once 
we know the probability distribution of the noise and that the constants C, 9 can be 
explicitly estimated from the hypotheses. 

We will assume throughout that the polynomial noise X is Gaussian and centered, the 
real-valued random processes 

Xi{-), . . . ,X„i{-) 

defined on R™ are independent, with covariancc functions 

i?^'(s, t) := E{X,{s)X,{t)) (^ = 1, . . . , m) 

depending only on the scalar product {s,t), that is: R'^'{s,t) — (5^''((s,t)), where 

di 

qW(u) = ^4'\'=, uem.{i = l,...,m). (6) 

k=0 

In this case, it is known that a necessary and sufficient condition for (5^*^((s,<)) to be 
a covariance is that c^ > for all k = 0,. . . ,di and the process Xi can be written as 

\\3\\<d. 

where the random variables a^* are centered Gaussian, independent and 

Var(af)=c«M: (* = 1, . . . ,m; ||j|l < d.) 

(for a proof, see, e.g., [1]). 

The Shub-Smale model (2) corresponds to the particular choice 

ci'^=(t) (fc = 0,l,...,d,) 
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)(0 



(u) = (1 + m)* (i = l,...,m). 



(7) 



We will use the following notations: 

Qu ,Quu denote the successive derivatives of Q''^. We assume that Q^^\u),Qu (u) do 
not vanish for u> 0. Put, for a; > 0: 



(0 



qiix) := 
nix) := 



Q_ 



hi(x) := 1 + a; 



qt{x)' 



(8a) 
(8b) 
(8c) 



In (8a) and (8b), the functions in the right-hand side arc computed at the point a;. 
In [1] the following statement was proved: 



Theorem 1. For any Borel set V C 



EiN-^iV)) 



iorel set V C K'" we 


tiave: 


1 ^(^\ 


I 


771 



1/2 



EdWtV^At 



(9) 



where 



En{x)^E\ 



E'^'^^)^' 



1/2n 



and ^1, . . . ,^m denote independent standard normal random variables. 

Remark. In fact. Theorem 1 is a special case of a general theorem (see [1]), in which the 
covariancc function of the random field is invariant under the action of the orthogonal 
group, and not only a function of the scalar product. 

Before the statement of our main result. Theorem 2 below, we need to introduce some 
additional notations and hypotheses. 

We will assume that each polynomial Q^*' docs not vanish for u>0, which amounts 
to saying that for each t the (one-dimensional) distribution oi Xi(t) does not degenerate. 
Also, Q*^*) has effective degree di, that is, 



^>0 



(i = l,...,m). 
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An elementary calculation then shows that for each polynomial Q*-*\ as u^ +00: 

^ 1 1 
h,{u) ~ -^ . (10b) 

Since we are interested in the large m asymptotic and the polynomials P, Q can vary 
with m, we will require somewhat more than relations (10a) and (10b), as specified in 
the following hypotheses: 

(Hi) hi is independent of i (i = 1, . . . ,to) (but may vary with m). We put h = hi. Of 
course, if the polynomials Q*^*) do not depend on i, this hypothesis is satisfied. 
But there are more general cases, such as covariances having the form (5(m)'' 
(i = l,...,m). 

(H2) There exist positive constants Di,Ei (i = 1, . . . , 771) and q such that 

0<Di-(l + u)qAu)< — and {l + u)qi{u)>q (11) 

1 + u - 

for all u > 0, and moreover 

max Di, max Ei 

l<i<m l<i<m 

are bounded by constants D,E, respectively, which arc independent of m; q is 
also independent of m. 

Also, there exist positive constants h, h such that 

h<{l + u)h{u)<h (12) 

for u>0. 

Notice that the auxiliary functions g^, r^, /i (i = 1, . . . , m) will also vary with m. To sim- 
plify somewhat the notation, we are dropping the parameter m in P, Q, qi, ri, h. However, 
in (H2) the constants h, h do not depend on m (see the examples after the statement of 
Theorem 2 below). 

With respect to (H2), it is clear that for each i, qi will satisfy (11) with the possible 
exception of the first inequality, and (1 + u)h[u) < h for some positive constant h, from 
the definitions (8a), (8c). (10a), (10b) and the conditions on the coefficients of Q*^*'. 
However, it is not self-evident from the definition (8c) that h{u) > for u > 0. This will 
become clear in the proof of Theorem 2 below. 

A second set of hypotheses on the system concerns the relation between the "signal" 
P and the "noise" X , which roughly speaking should neither be too small nor too big. 

Let P be a polynomial in m real variables with real cocfhcicnts having degree d and 
Q a polynomial in one variable with non-negative coefficients, also having degree d, 
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Q{u) = J2k=o ^kU^ ■ We assume that Q does not vanish on u > and Cd > 0. Define 



H{P,Q):^ sup (l + ||t||) 



tei 



V 



P 



^/QW¥) 



it) 



K{P,Q):= sup Ul + ||if) 



d / P 

d~p\7QWF) 



it) 



where ^ denotes the derivative in the direction defined by -Aw, at each point t 7^ 0. 
For r > 0, put: 



L{P,Q,r):= inf 



Pit)' 



t\\>rQi\\tP)- 

One can check by means of elementary computations that for each pair P, Q as above, 
one has 

H{P,Q)< 00, K{P,Q)< 00. 

With these notations, we introduce the following hypotheses on the systems P,Q; as m 
grows: 

(H3) 

=1 



Ar^ 



1 ™ 

-•E- 



= 0(1) asm ^+00, 



B„ 



1 ^ i^^(P.,QW) 
— ■ > ^ ~ 0(1) as ?Ti ^ +00. 



(13a) 
(13b) 



(H4) There exist positive constants ro,£ such that if r > vq: 

L{P,,Q'-'\r)>e foraUi = l,.. 
Theorem 2. Under the hypotheses (Hi), . . . , (H4), one has 

E{nP+^)<C9"'E{N^), 
where C, 9 are positive constants, < 9 <1. 



(14) 



1.1. Remarks on the statement of Theorem 2 

1. In fact, wc will sec in the proof of the theorem how one can get explicitly from the 
hypotheses first the value of 9 and then mo and the constant C in such a way that 
whenever m > ttiq, inequality (14) holds true. 
A possible choice is as follows: 
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• Choose ro from (H4), 

''0 -f/2] n 1 + 1 



• Let us put Fi = Ei/Di (z = 1, . . . , m) and F = max{Fi, . . . , Fm}- 
From the hypotheses, one has F < E/q. Let t > such that: 

F 1 1 

13 < 97X72- i^^i 



1 + rV^ 2 1 + r^ 
Choose mo (using (H3)) so that \i m> ttlq one has: 

g[mA„/g+rnS„/(hg)]/2^ni^ ^ 
^2-,2 \ (m-l)/2 



(16) 



"'.i+Xgj <V7^- 



Then, (14) is satisfied for m > mo, with 



/lyr 



-r, 



2 



C = 30--^^ — — ^. (17) 

h ra 

2. It is obvious that our problem does not depend on the order in which the equations 

P^{t)+X,{t)^0 {i^l,...,m) 

appear. However, conditions (13a) and (13b) in hypothesis (H3) do depend on the 
order. One can state them by saying that there exists an order z = 1, . . . ,m on the 
equations, such that (13a) and (13b) hold true. 

3. Condition (H3) can be interpreted as a bound on the quotient signal over noise. 
In fact, it concerns the gradient of this quotient. In (13b) the radial derivative 
appears, which happens to decrease faster as ||tj| -^ 00 than the other components 
of the gradient. 

Clearly, \iH{P„ Q'*)), K{Pi, Q'*)) are bounded by fixed constants, (13a) and (13b) 
are verified. Also, some of them may grow as vi -^ +00 provided (13a) and (13b) 
remain satisfied. 

4. Hypothesis (H4) goes ~ in some sense - in the opposite direction: For large values 
of ||i|| we need a lower bound of the relation signal over noise. 

5. A result of the type of Theorem 2 can not be obtained without putting some re- 
strictions on the relation signal over noise. In fact, consider the system 

P^{t)+(jX,{t)^0 (z = l,...,m), (18) 

where cr is a positive real parameter. For generic P, as cr | the expected value of 
the number of roots of (18) tends to the number of roots of Pi(i) = (i = 1, . . . , m). 
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In this case, the relation signal over noise tends to infinity. On the other hand, if we 
let (T —> +00, the relation signal over noise tends to zero and the expected number 
of roots will tend to E{N^). 

2. Some examples 

2.1. Shub Smale 

In the Shub-Smalc model for the noise, Q*-*^ is given by (7). Then, 

qi (u) = -— -^ , hi (u) ^ h{u) 



1+u " ' " ' 1+u 

We assume that the degrees di arc uniformly bounded. So, (Hi) and (II2) hold true. 
Of course, conditions (H3) and (II4) also depend on the signal. 
We are going to give two simple examples. Let 

where di is even and r is positive and remains bounded as m varies. One has: 
d f P, V, c;jtr--i+d,r'^-|l^ll ^ d,{l + r^^) 

/ p, \ d.wtw^^^'^ + d^' 

ddl+r'i^) 



which implies 



'(JwK" 



< 



(l + ||t|12)3/2- 



Again, since the degrees di, . . . , dm arc bounded by a constant that does not depend on 
m, (H3) follows. (II4) also holds under the same hypothesis. 

To illustrate a numerical example of the values of 6, C and m^, let us assume that the 
radius r = 1. Then: 

• We can choose tq = 2. It turns out that t = 9/25 if Z) < 4 and ^ = (2^ - 1)^5^ if 
D>b. 

• 6'=(3 + 2V2)/6. 

• C = 15 v5 and toq can be chosen so that \l m> mo, then 

gCi^ci + 1/2 < ^m ^^^^^ ^ g^2 ^ ^ — + - > 1. 

4V2 2 
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Notice that an interest in this choice of the P^'s hes in the fact that obviously the 
system Pi{t) = {i = 1, . . .,m) has infinite roots (all points in the sphere of radius r 
centered at the origin are solutions) , but the expected number of roots of the perturbed 
system is geometrically smaller than the Shub-Smale expectation, when m is large. 

Our second example is the following: Let T be a polynomial of degree d in one variable 
that has d distinct real roots. Define: 

P,(ii,...,t„i) ^T{ti) (i = l,...,m). 

One can easily check that the system verifies our hypotheses, so that there exist C, 9 
positive constants, < < 1 such that 



e{np+^) < ce^'d" 



72 



where we have used the Shub-Smalc formula when the degrees are all the same. On the 
other hand, it is clear that N^ ~ d™ so that the diminishing effect of the noise on the 
number of roots can be observed. A number of variations of these examples for P can be 
constructed, but wc will not pursue the subject here. 

2.2. Q^^^ = Q, only real roots 

Assume all the Q^*-' arc equal, Q^*' = Q and Q has only real roots. Since Q docs not 
vanish on u > 0, all the roots should be strictly negative, say 

-«!,. ..,-ad, 

where 

< ai < a2 < • • • < ctd- 

With no loss of generality, we may assume that ai > 1. If this were not the case, we 
perform a homothecy of the space M™, centered at the origin and with factor equal to 
y/ai- The number of roots remains unchanged, and the new Q has ai > 1. 

We will assume again that the degree d of Q is bounded by a fixed constant d (one 
should take into account that Q may vary with m), as well as the roots 

ak<a (fc = l,...,d) 

for some constant a. A direct computation gives: 

d 

ak 



q^{u) = q{u)^y2^, , hiiu) = h{u) ^ ——Y] 

^ — ' 7/ 4- rvi_ rill \ ^ — ^ 



k=\ 



M + afe' g(li) ^ (u + afc)2 



One verifies (11), choosing Di = d, Ei = dia^ — 1). Similarly, a direct computation gives 
(12). 
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Again let us consider the particular example of signals: 

where di is even and for each i — 1, ... ,771 and r is positive and remains bounded as m 
varies. 






<d,(S + r'*0 ^ 



(l + ||t||2)3/2 



so that K{Pi,Q'-'^') is uniformly bounded. A similar computation shows that H{Pi,Q^^'') 
is uniformly bounded. Finally, it is obvious that 

L{P,,Q^'\r)> 



1 + a 



for i ~ I, . . . ,m and any r > 1. So the conclusion of Theorem 2 can be applied. 

One can similarly check that the second polynomial system in the previous example 
also works with respect to this noise. 

2.3. More general examples 

Assume that the noise has covariance with the form 

Q^'\u,v,w)^[Q{u)]'' (z = l,...,m), 

where Q is a polynomial in one variable having degree v with positive coefficients, Q{u) = 
Sfe=o bkU ■ Q may depend on m, as well as the exponents li, . . . ,1^. Notice that di — v-li 
{i = l,...,m). 

The Shub-Smale case corresponds to the simple choice Q(u) = 1 + u, li — di (z = 
l,...,m). 

One has: 

Q{u) 

Q'Hu)-Q{u)Q"{u) 



hi{u) = h(u) = 1 — u 



Q{u)Q'{u) 



so that (Hi) is satisfied. 

We will require the coefficients bo, . . . ,bi, of the polynomial Q to verify the conditions 

bk< ''~ j^^ bk^i (fc = l,2,...,i/). 
Moreover, we assume that 
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are bounded by a constant independent of m and there exist positive constant b, b such 
that 

b<bo,bi,...,b^<b. 

Under these conditions, one can check that (H2) holds true, with Di = di (i = 1, . . . ,77i). 
For the relation signal over noise, conditions are similar to the previous example. 
Notice that already \i v ~2 and we choose for Q the fixed polynomial: 

Q{u) = l + 2au + bu^ 

with < a < 1, v6 > a > 6 > 0, then the conditions in this example are satisfied, but the 
polynomial Q (hence Q"^') does not have real roots, so that it is not included in Example 
2.2. 



3. Proof of Theorem 2 

Proof of Theorem 2. Let 

P,{t)+X,{t) 



ZJt)^^^ '' ' (7 = 1,. 



and 




Clearly, 

for any subset V of W^ . 

Clearly, the Gaussian random fields {Zj{t):tE M™} (j = l,...,m) are independent 
and 

V^T{Z]{t))=E(Z]it)) = l (19) 

for all j = 1, . . . , m and aU t e M™, where Zj{t) = Zj{t) - E{Zj{t)). 

Differentiating in (19) with respect to ta{a — 1, ■■■ ,m) we obtain that E{^^Zj{t)Zj{t)) 
0, (j = 1, . . . , m). Since the joint distribution is Gaussian, this implies that \7 Zj{t) is in- 
dependent of Zj{t), that is, Z'{t) and Z{t) are independent. 

We apply the Rice formula (see [2] for a complete proof) to compute E{N^{V)), that 
is: 

E{N^iV))= f E{\dct{Z'm\Z{t)^0)-pz^t^iO)dt, 
Jv 
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where pj(-) denotes the density of the probabihty distribution of the random vector ^, 
whenever it exists and f'{t) : R" -^ M'" the derivative of the function / : M™ ^ R" at the 
point t. Because of the independence between Z'{t) and Z{t), we can erase the condition 
in the conditional expectation, obtaining: 



E{N^{V))^ / E{\dct{Z\tm 



y (27T)"/2 



(20) 



Our next problem is the evaluation of E{\ dct(Z'(f))|). 
A direct computation of covarianccs gives: 



fdZ dZ, 



92 



^^))-^-a;^^''(^'*)l- 



2^J- + I „ /IUI|2^ 



5^J[n{\\t\\)to.tf3+q^{\\t\\)5c|}\ 



for z, j, a,/3 = 1, . . . ,771, where the functions qi^Ti have been defined in (8a) and (8b). 

For each t ^ 0, let Ut be an orthogonal transformation of R™ that takes the first 
element of the canonical basis into the unit vector -An ■ Then 



Varf i^Si V Diag(MI|i|r)a, 
where we denote the gradient VZj(i) as a column vector. 



.,1), 



(21) 



^f">^ 



VZ,{t) 



dZ 



Kmt^'V 



Diag(Ai, . . . , Am) denotes the mxm diagonal matrix with elements Ai,. 
diagonal and the function h has been defined in (8c). So we can write 



, Am in the 



UtVZ,(t) 



Cj+ttj (j = l,...,777). 



where Q is a Gaussian centered random vector in R™ having covariance given by (21), 
Ci : ■ • ■ I Cm sie independent and aj is the non-random vector 



UtyiP,{t)/VQ^^H\\m) 



( 



^li 



(,? = i,. 



(22) 
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Denote T as the m x m random matrix having columns r]j = Ci + Q^i (j = I7 • ■ • i 't^)- We 
have 

m 

\det{Z'm^\dct{T)\-l[{q,{\\tr)r^' 
so that 

m 

E{\ det(Z'(t))|) = £;(! det(r)l) . n(*(lltf ))'/'■ (23) 

The components C,ij (i = 1, . . . , m) of Cj are Gaussian centered independent and 
Var((^ij ) = 1 for i = 2, . . . , m; j = 1, . . . , m 

Var(Ci,) = MIlif ) forj = l,...,m. 

Put 



so that 






\dctiT)\^^m¥)-\dctif)\, 



where T is the random matrix having columns 77^ = tjjj +aj with ipi, . . . , ■(/;,„ i.i.d. standard 
normal in M™. 

I det(r)| is the volume of the parallelotope generated by 771, . . . ,7),„, that is, the set of 
vectors in R™ which can be expressed as 

m 

with < Ci < 1 (1 = 1,..., m). Hence, 

|det(T)| = y^Om • ll^™ll • n div,,S,), (24) 

where Sj {j = 1, . . .,m—l) stands for the subspace of M™ generated by fjj^i ,.. .,fjm and 
d is Euclidean distance. Notice that for j = 1,...,to— 1, using invariance of the standard 
Gaussian distribution under the orthogonal group, 

E{d{fjj,S,) I %+i, . . . , ,7™) = E{U, + d,\\j), (25) 
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where || • \\j denotes Euclidean norm in R-' (|| • || = || • ||m), Ci is standard normal in W 
and dj is the orthogonal projection of aj onto the orthogonal complement of Sj in R™ 
(which is identified here with W), so that ||aj||j < ||i3j||- Introduce now the function 

where c€W is non-random. It is clear that jj is a function only of \\c\\j, which is in fact 
increasing and 

7,(c)<(l + i|c|||l)7,(0). 

(See the auxiliary Lemma 1 after this proof.) So, if we denote by Cj a non-random vector 
in W such that \\cj\\j = \\aj\\j, it follows from (24) and (25), by successive conditioning 
that 



i?(idetri)<v^^MF)-n^(ii^^+^^ii^)^VMim- n^^-(o) n(i 



i=i 



Vi=i 



Vi=i 



,1_ 

'2j 



Using (20) and (23) we get: 



E{N^) < 



(2nynl2- 



t{v^om-(n «■<««) 



1/2 



X exp 






(26) 



dt, 



where 



Our final task is to obtain an adequate bound for the integral in (26). For 7 = 1, . . . ,77i 
(use (H2)): 



>^ijl 



^jm¥wW) 



d mtW) 



^pVWWF) 



< 



1 



-.K{P„Q 



Uh 



and 



V^MW) "VI 



Then, if we bound W&jW^ by: 



|fi,||'<|5i,f + ||a,||2 
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we obtain 

which implies (using (H3)) 

Y] IICjIl ■ • - < -rnAm + —772Bm. 

Replacing in (26) we get the bound: 

E(N^) < SmH,n, (27) 

where 

"■ \ ,g(mA„/g+mS„/(?ig))/2 ^28) 





1 An 




r 




^7^(m+l)/2 ^2 



„ / m \ 1/2 

/ (ll'^^di^ii')) v^^MF)i?(!i6nii)e-'i:".^'(*)^/'3-(ii*in)/.di. 



(29) 



The integrand in (29) is the same as in formula (9) giving the expectation in the centered 
case, except for the exponential, which will help for large values of ||<||. 
Let us write Hm as 

where Hm (j) corresponds to integrating on ||t|| < r and Hm (r) on ||t|| > r instead of 
the whole R™ in formula (29). We first choose r — tq so that, using H4: 

H^n\ro) < e-' ™/2^(A^-^). (30) 

We now turn to Hm (r). We have, bounding the exponential in the integrand by 1 and 
using hypothesis (H2 ) : 

where am-i is the (to — l)-dimensional area measure of S""^^. The integral in the right- 
hand side is bounded by 

/ 2 \ (m-l)/2 

2\l + r^ 
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Again using (H2) and formula (9), we have the lower bound: 



r(W2) ,1/2 



+CX3 



nu 



D, 



E, 



2\2 



\\l+P^ (l+p2) 



1/2' 



(l+p2)l/ 



TH'^P 



y27t('"+i)/' 



^/i^/^i?(ii6„ii) nA^'M-"-! 



vj=i 



(l+p2)('»+ 



4 = 1 ^ 



F, 



1/2 



^7g^i'"^«-j'(n«-"^)--. 



dp 



rl + l 



+ CK3 



1)71^^' 



(l+/?2)(™+l)/ 

where r has been chosen to satisfy (15). 

To get a lower bound for the last integral, a direct integration by parts shows that: 

+00 m~l ^-2 

(l+p2)(m+l)/2^/'>^' 

which implies 



+00 



e-2 nf T-'rl \ ("-1^2 



/rro (l+p2)(™+l)/2'^^> V^ 2 (^ 1 + T2r2 

Then, choosing toq as required by condition (16), we get for m > ttiq: 



m^Kro)<C^ 



rl + 1/2 



m/2 



E{N^) 



(32) 



where 



Ci = 7tc^ - 



/.V/^v/rT72 



^"0 



For the remainder, we must put together (27), (28), (30) and (32). One easily checks 
now that with the constants given by (17), the inequality of the statement holds true. D 



3.1. Auxiliary lemma 

Lemma 1. Let 7 : M'' ^ M, fc > 1 be defined as 

7(c) = i?(iie+c|i), 
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where S^ is a standard normal random vector in Mr, and c E R'^ (\\ ■ \\ is the Euclidean 
norm in R'^j. Then 

(i)7(o) = V2^:f±^. 

(ii) 7 is a function of \\c\\ and verifies: 

7(c) < 7(0) (l + ^llc|p). (33) 

Proof, (i) In polar coordinates 7(0) = (p)jy2 Jn °° p^c~'' 1'^ . Using the change of variable 
u = p^ , we obtain the result. 

(ii) That 7 is a function of ||cj| is a consequence of the invariance of the distribution 
of S, under the isometrics of R*"'. For fc = 1, (33) follows from the exact computation 

7(c) = y27^e-^/2 + c|'^^e--^/^ 

and a Taylor expansion at c = 0, which gives 

1 2 



7(c)<v/2/7t(1 + -c 

For k>2, we write 

7(c) = £;([(ei+a)2+e| + --- + 4']'/')-G(a), 

where a = \\c\\ and ^1, • ■ • ,^fc are independent standard normal variables. Differentiating 
under the expectation sign, we get: 

G'ia)=E[ ^^+" 



so that G"(0) = due to the symmetry of the distribution of ^. 
One can differentiate formally once more, obtaining: 

For the validity of equality (34) for fc > 3 one can use that if d > 2, -n^ is integrable in 

R'* with respect to the Gaussian standard measure. For k = 2 one must be more careful 
but it holds true. The other ingredient of the proof is that one can verify that G" has a 
maximum at a = 0. Hence, on applying Taylor's formula, we get 

G{a)<G{0) + ^a^G"{0). 
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V2 T{{k+l)/ 
k r(fe/2) 



Check that G"(0) = ^^^^757^ which, together with (i), gives: 



C?"(0) _ 1 
G(0) " fc' 

which iinphes (ii). D 
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